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Abstract 

The group-theoretical method used to derive horizontal symmetry from neutrino mixing is re- 
viewed and expanded. Some misunderstanding in the literature regarding the result is clarified. 
The method used previously to find vacuum alignments of £4 is applied to compute those of A4 
and 53. A study of effective theories based on these three groups shows that in each case there 
are just enough free parameters to fit all the masses and the remaining mixing parameters. This 
places constraint on dynamical models because effective theories are just dynamical models with 
the right-handed fermions integrated out. How quarks may fit into this scheme is briefly discussed. 
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I. INTRODUCTION 



The usual method to model a mass matrix proceeds from the top down. First, a horizontal 
symmetry group Q is picked, an irreducible representation (IR) is assigned to each left- 
handed (LH) and each right-handed (RH) fermion in the theory, and a number of Higgs-like 
scalar bosons is created to couple invariantly to the fermions. Then, Higgs expectation 
values are introduced to break the symmetry, from which mass matrices are computed. The 
mass matrices so obtained depend on the Yukawa couplings, the mass parameters, and the 
Higgs expectation values. These parameters are chosen to fit the observed fermion masses, 
the mixing angles, and the phases. Such a construction is often applied to the leptons alone, 
but sometimes it is also carried out in grand unified theories containing both leptons and 
quarks. In the latter case frequently only approximate fits are obtained. 

A model containing an equal or a larger number of tunable parameters than available 
experimental quantities can be expected to fit all the data, even without a horizontal sym- 
metry. Thus the introduction of a horizontal symmetry is justified only when it leads to 
unambiguous predictions to allow a smaller number of parameters to be used to fit the data. 
Many existing models do have fewer parameters, some much fewer, but if they fit the ex- 
perimental quantities only approximately, then it is unclear what and where exactly are the 
savings coming from the symmetry assumption. 

In this article we review and expand on a complementary approach [1], using only left- 
handed fermions and start from the bottom up. In this approach, symmetry is directly 
connected to mixing, so the gain from symmetry is clear. The simplification of using only 
LH fermions comes about because low energy data such as masses and mixing parameters 
can all be determined from the effective mass matrices connecting only LH fermions. 

The simplification also brings about its limitation, in that without the RH fermions 
local dynamical theories cannot be constructed. In other words, although the connection 
between mixing and the LH effective theory is straight forward and unique, the next step 
of constructing a dynamical model is not. Nevertheless, this analysis tells us what would 
happen to every dynamical model after its RH components are integrated out. To the extent 
that low energy data can all be determined from the LH effective theories, they already carry 
most of the information reachable by present-day experiments. 

Sec. 2 contains a review of the general method and the previous results. This method 
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is used to compute the vacuum alignments of A 4 and S3 in Sec. 3. The results are listed 
in Table 1, together with the alignments of £4 previously obtained and the IR from which 
these vacuum alignments are calculated. 

A detailed discussion of the LH effective theories for these groups is carried out in Sec. 4. 
The number of independent parameters is determined in each case, with the result tabulated 
in Table 2. It turns out that in all the cases considered, there are always just enough 
parameters to fit all the masses and the undetermined mixing parameters. 

In principle, this bottom-up group-theoretical approach can be applied to the quark sector 
as well, but the symmetry so obtained is very different from the symmetry in the lepton 
sector, because their mixing matrices are vastly different. For reasons already discussed [1], 
we prefer to determine the horizontal symmetry from the lepton sector. The question of 
how quarks can fit into this scheme will be briefly discussed in Sec. 5. 

The remainder of this section is devoted to addressing some confusions found in the 
literature regarding our previous results [1], summarized below. 

1. If the neutrino mixing matrix is tri-bimaximal (TBM) [2], then the minimal left-handed 
(LH) symmetry group 1 in the lepton sector is S4; 'minimal' means that every other 
LH symmetry group contains S4 as a subgroup. 

2. If only trimaximal mixing taken from the second column of the TBM is assumed in 
the mixing, then the smallest LH symmetry group is A4. It is not known whether this 
group is minimal or not. 

3. If only bimaximal mixing taken from the third column of the TBM is assumed, then 
the smallest LH symmetry group is S3. It is not known whether this group is minimal 
or not. 

4. If only the first column of the TBM is assumed, then the minimal symmetry group 
is again £4. In principle the mixing matrix so obtained need not be TBM, because 

1 This is the horizontal symmetry group for an effective theory in which only LH fermions and composite 
Higgs appear. Such an effective theory is obtained from a dynamical model by integrating out its RH 
fermions. The resulting effective Hamiltonian is generally non-local, but that does not matter because we 
are only interested in its horizontal-symmetry structure. This group was referred to as the natural sym- 
metry group in [1], but since there is objection to that name [4], I will use a more descriptive terminology 
here and call it the LH symmetry group. 
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neither trimaximal mixing nor bimaximal mixing is explicitly assumed. 

Starting from a symmetry group Q so derived, one can compute the Higgs alignment 
needed to regain the original mixing. This has been carried out in S4 [1], and it will also be 
done for A 4 and £3 in Sec. 3. 

Obviously additional soft breaking or some other vacuum alignment can be imposed on 
Q so that the original mixing input is not recovered. Conversely, one may fit all the experi- 
mental data including mixing without even any symmetry as long as a sufficient number of 
parameters are provided. The present approach simply tells us that certain predictions on 
mixing can be obtained if an appropriate LH group Q and the appropriate vacuum align- 
ments are used, and that this correspondence is one to one. 

If Q is the LH symmetry group, one can always construct a local dynamical model with 
Q as the dynamical symmetry, giving rise to the same prediction on mixing, but there is no 
claim to the converse. In other words, it is conceivable that a dynamical model with some 
other symmetry group can lead to a LH effective theory with an enhanced LH symmetry 
group Q. With these remarks, one finds that the two 'counter examples' in [3] really do not 
contradict the results cited in 1 to 4 above. 

II. THE CONNECTION BETWEEN MIXING AND THE HORIZONTAL GROUP 

In this section we review the relation between the mixing matrix and the minimal LH 
symmetry group Q [1]. The relation comes from the recognition that Q is the minimal group 
generated by a set of residual symmetry operators that can be computed from the mixing 
matrix. Conversely, given Q, the mixing matrix can be retrieved if certain computable 
vacuum alignments (direction of Higgs expectation values) are used. The rest of this section 
contains a sketch of the details to establish terminologies and notations. It also explains 
how the familiar groups S4, At, and S 3 for leptons arise in this context. 

Let Mi (i = u, d, e, v) be the Dirac mass matrices for the up-quarks, down-quarks, 
charged-leptons, and neutrinos respectively, connecting LH to RH fermions, and Mjy be 
the Majorana mass matrix for the heavy ('right-handed') neutrinos. Since we only deal with 
LH fermions, the mass matrices we need to consider are the LH to LH effective mass matrices 
Mi, defined by M t = M { M\ (i = u, d, e) and M v = M r M x : M', . Let Ui (i = u, d, e, v) be the 
unitary matrices that render UjMiUi (i = w, d, e) and UlM v U v diagonal. Then the mixing 
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matrix is Uckm = U^Ua for the quarks and Upmns = F\U V for the leptons. In particular, 
in the basis where Md and M e are diagonal, which we shall adopt from now on, Uq KM = U u 
and Upmns — U v . In most of the following we deal explicitly with the lepton sector, but 
the quark sector can be treated analogously. 

Let F and G be respectively unitary symmetry operators for charged leptons and neu- 
trinos, meaning that F^M e F = M e and G T M U G = M v . Then it follows that F and M e 
have common eigenvectors, given by the columns of U e — 1, and G and M v have common 
(pseudo-)eigenvectors, given by the columns of U u = Upmns- Being unitary, the eigenvalues 
of both F and G are unimodular, but on account of the Majorana character of the neu- 
trinos, the eigenvalues of G must also be real, and by convention one of them is chosen to 
be +1 and two of them chosen to be —1. With this convention, three symmetry operators 
G a (a = 1,2,3) can be constructed for each mixing matrix U v , in which the eigenvector 
with +1 eigenvalue is taken from the ath column of U v . These operators satisfy G 2 a = 1 and 
G a Gb = G c if a, b, c are different, so only two of the three are independent. In contrast, F 
remains uncertain because it may have arbitrary eigenvalues, as long as they are unimodular. 

Abstractly, the group generated by G a is always Z 2 x Z 2) whatever the mixing matrix is, 
and the group generated by F is always a subgroup of U(l). Obvious that information is 
useless [1, 3] if we do not know how this subgroup in the charged- lepton sector is intertwined 
with the group Z 2 x Z 2 in the neutrino sector. To be useful we need the detailed form of F 
and G a determined from the mixing matrix to establish the horizontal LH symmetry group 

Q. 

The symmetry operators F and G shall be referred to as residual symmetry operators, 
and we will use the notation 1Z = (F, G) to denote them. We know that F ^ G or else there 
would be no mixing. The minimal horizontal LH symmetry group Q is then the minimal 
group generated by the operators F and G: Q = {F, G}. 

In principle there could be many residual operators of type F, but to find minimal groups 
we may confine to only one of them. If we demand Q to be a finite group, we must have 
F n = 1 for some finite n. If we want to be able to reconstruct the mixing matrix from the 
residual symmetry operators, we must require F to be non-degenerate, in the sense that 
no two eigenvalues are identical, hence n > 3. Even with these restrictions we still have 
an infinite choice of F, by picking any n > 3, any three distinct nth root of unity as its 
eigenvalues, and one of the 3! = 6 permuted arrangements of eigenvalues in the matrix F. 
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Correspondingly, there may be many minimal LH symmetry groups Q. 

We shall return to the choice of F after discussing the experimental implications on G a . 
In the neutrino sector, we shall first assume Upmns — U v to be of the tri-bimaximal (TBM) 
form 

/ 2 V2 \ 



V6 



(1) 



-1 V2 
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which when compared with experimental data is accurate to within one standard deviation. 
Using the recipe explained above, we can compute from that the three residual symmetry 
operators to be 
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Now return to the choice of F. For n — 3, the three distinct eigenvalues are the cubic 
roots of identity, l,u — e 27 ™/ 3 , and to 2 . There are six distinct F's corresponding to the 
3! permutations, but only two of them generate distinct groups Q. The one generated by 
F\ = diag(l, cu, uj 2 ) is Q = {Fi, G2, C3} = 5*4, the permutation group of four objects, and 
the one generated by F 2 = diag(u;, 1, uj 2 ) is Q = {F 2 , G 2 , G 3 } = 3^4, consisting of Sa^Sa, 
and cu 2 ^. It contains S4 as a subgroup. One representation of F ± , G±, G 2 , G s in ^4 is 
F 1 = (123), G 1 = (14), G 2 = (14)(23), G 3 = (23). This representation will be useful in 
seeing how subgroups of S4 are generated. 

For any choice of F with n > 3, it can be shown that if Q = {F, G 2 , G3} is a finite 
group, then it always contains S4 as a subgroup. Hence 64 is the minimal symmetry group 
for TBM. For example, PSL 2 (7) = £(168) is one of these groups [5] containing S4 as a 
subgroup. 

The tri-bimaximal matrix implies a vanishing reactor angle #13, with an ineffective CP 
phase. Since the experimental error of this angle is still large, it is reasonable not to be so 
committed at the present, which means that the bimaximal character appeared in the third 
column of (1) should be abandoned. We may still assume the second column (trimaximal 
mixing) of (1) to be valid [6], in which case the first and third columns of U u are parameterized 
by two real parameters. With that, the only known symmetry operator in the neutrino 
sector is G 2) and the horizontal group generated by {F 1: G 2 } can be seen to be Q = A 4 , 
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the subgroup of S4 with even permutations. If the two parameters in the mixing matrix are 
adjusted so that #13 = 0, then an accidental symmetry exists to enlarge the symmetry from 
A 4 to S4. The LH symmetry group generated by other F's and G 2 are not presently known. 

Similarly, if we retain the first column of (1), then Q = {Fi,Gi} = S4, but in this case 
the reactor angle does not necessarily vanish. This illustrates the fact the symmetry alone 
cannot determine the mixing matrix. What is needed in addition is the alignment of the 
composite Higgs, to be discussed in the next section. 

If for some reason we only want to retain the third column of (1), namely, the bimaximal 
mixing, then it is well known that the reactor angle vanishes and the atmospheric angle 
is maximal. In that case the solar angle is a parameter, and the symmetry group is Q = 
{F 1 ,G 3 } = S 3 . 

III. VACUUM ALIGNMENT 

Given a LH symmetry group Q, the original mixing can be recovered by keeping the 
residual symmetry 1Z intact in the spontaneous breaking of Q. This in turn can be achieved 
by having the correct vacuum alignments computed as follows. 

We start with an effective Hamiltonian H, constructed from LH fermions and composite 
Higgs fields. It can be thought of as the result of a dynamical Hamiltonian H with the RH 
fermions integrated out. The Higgs field in H are often bilinear in the Higgs fields in H, 
hence the name 'composite'. H may not be local, but it is invariant under the LH symmetry 
group Q. We retain in H only the Yukawa terms, because only those contribute to the mass 
matrices. Mass terms may be present but they are special cases of Yukawa terms coupled 
to a singlet Higgs possessing non-zero expectation values, so they will not be considered 
separately. The composite Higgs coupled to the LH charged leptons will be denoted by $, 
and those coupled to the LH neutrinos will be denoted by See Sec. 4 for details. 

The condition for 1Z = (F, G) to remain a symmetry of H after the Higgs acquire their 
expectation values is 

for every residual symmetry operator G and for every IR (A) of the group Q. F^ and G^ 
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are the IRs of F and G. 

For Q = {Fi, G 2 , G3} = 6*4, these alignments have been computed, and the result is listed 
in Table 1. The other three cases in Table 1 are computed similarly using (3) and the IRs 
of F^ and G^ A \ also listed in Table 1. Note that the 3-dimensional representation of S3 
is always reducible. In that table, Q = {Fi,G±} = £4 is denoted £4, to distinguish it from 
Q = {Fi, G 2 , G3} which is denoted as £4. It has a different residual symmetry 71, so its 
Higgs alignments are not the same as £4. 

In Table 1, the IRs are listed in the second row, with numerals indicating the dimension 
of the representation, and — indicating the absence of a particular IR in that group. For 
easiness of printing, column vectors are written as row vectors in the table, and (1, 1, 
is any vector orthogonal to (1, 1, 1). The overall normalization of each vacuum alignment 
is arbitrary because (3) is linear, 'dr' in the first row stands for 'defining representation', 
which is the representation the LH fermions belong to. I apologize for the convention 
adopted here in which the defining representation of £4, for which det(Gj) = +1, is called 
3' rather than 3. The other symbols are, u = e 2 ™/ 3 , <7\ is the Pauli matrix [[0,1], [1,0]], 
and g a are the G a defined in eq. (2), namely, g\ = [[1, — 2, — 2], [— 2, — 2, 1], [— 2, 1, — 2]]/3, 
g 2 = -[[1, -2, -2], [-2, 1, -2], [-2, -2, l]]/3, and g 3 = -[[1, 0, 0], [0, 0, 1], [0, 1, 0]]. Note that 
in the 3' IR is (0,0,0) for £4 but (1, 1, l)j_ for £4. It is this difference that in principle 
allows £4 to accommodate non-zero 9 13 , though the discussion in the next section shows 
that this turns out to be impossible. 
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Table 1. Vacuum alignments ($) and of the horizontal symmetry groups At, 6*3, S4 
for leptons. See the text for an explanation of the symbols. 

IV. LOW ENERGY CONTENT 

As mentioned in the Introduction, the low-energy content (masses, mixing angles, and 
phases) of any dynamical theory is already contained in the LH effective mass matrices 
Mi (i = u,d,e,is), which can be read off from the LH effective Hamiltonian H. In this 
section we shall discuss the structure of H for every symmetry group Q discussed in the 
previous section. 
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The effective Hamiltonian can be written symbolically as 

H = J2 (*AeW® A + ^I^ A ) , (4) 

A 

with all couplings understood to be (/-invariant. The LH charged-lepton fields are denoted by 
eL, and the LH neutrino fields are denoted by vl. They belong to the defining representations 
r = 3', 3, (1, 2), 3', respectively for the groups Q = S 4 , A 4 , S 3 , S 4 (see Table 1). Only IRs A 
contained in the Clebsch-Gordan (CG) series r x r appear in the sum in (4). The effective 
mass matrices M e and M v can be read off from (4) using the Higgs alignments listed in 
Table 1. Since M e = MJ and M v = Mj, IR A that leads to an anti-symmetric contribution 
to M v should be dropped. Moreover, the parameters A^ have to be real though /ia m &y be 
complex. 

The construction of (4) is schematically shown in Table 2. The first and second rows 

together specify the defining representation r of the various groups, and the remaining rows 

indicate how the CG series r x r pans out. — denotes an IR that is not present in that group, 

and x denotes an IR that does not appear in the CG series r x r. A '0' indicates a term 

in which the expectation value ($) or vanishes so the corresponding coupling constant 

does not appear as a parameter. Each entry in S3 is shown as a 2 x 2 matrix, indicating the 

/lxl lx2\ 

Clebsch-Gordan series for 

\2xl 2x2/ 

The directions of the Higgs expectation values are given in Table 1, and their magnitudes 
can be absorbed into the Yukawa coupling constants A^ or /i^, hence the number of free 
parameters is just the number of Yukawa coupling constants, each of which is represented 
by a black dot, with the total number listed in the last column. An open circle indicates 
an antisymmetric contribution to M u that has been dropped. In the case of S3, the open 
circle actually represents the antisymmetric combination of 1 x 2 and 2x1. As mentioned 
above, the parameters in the <£> rows have to be real but the parameters in the \1/ rows may 
be complex. 
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Table 2. Coupling scheme for H in (4), with the total number of independent free parameters 
listed in the last column. See the text before the table for an explanation of the symbols. 

We want to find out whether mass matrices so obtained contain less, just enough, or more 
parameters than necessary to fit the low-energy data. In the charged-lepton sector, every <3> 
row of Table 2 contains exactly three real parameters to fill the three entries of the diagonal 
mass matrix M e , which is just enough to fit the three charged-lepton masses no matter what 

g is. 

The neutrino sector is more complicated. The three neutrino masses combining with the 
two Majorana phases form three complex masses, albeit with an arbitrary overall phase. 
As to mixing, since £4 automatically gives rise to the tri-bimaximal matrix, no additional 
parameter is required to fix the mixing, so the total number of complex parameters needed 
to fit the low energy data is 3, exactly what the effective Hamiltonian provides. In the other 
three cases, where only one column of the tri-bimaximal matrix is automatically obtained, 
we need one additional complex parameter to fully specify the mixing matrix [7]. Thus in 
those three cases the total number of complex parameters needed to fit the low energy data 
is 4. From Table 2, we see that is precisely what we have for A 4 and S3, but we are one 
short in S4. As explained below, it turns out in this case that #13 = automatically, so all 
the mixings are determined, and 3 is just the right number of parameters to fit the complex 
neutrino masses. 
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(5) 



Let us now look at the resulting neutrino mass matrices to understand what really hap- 
pens. For At, explicit calculation shows that 

/ fii + 2fi 3 fii" - fi 3 fir - fi 3 \ 
My = fii" - fi 3 fir + 2/i 3 fii - fiz 

\ fi r - fi 3 fii- fi 3 fir, + 2^3 J 

where fii = /ii/v^, fir = Hi'/V3, fii" = fii"/\^3, fi 3 = n 3 j\f§. This matrix is magic 
[8], whose row and column sums all equal to fii + fir + fii", hence the mixing matrix 
computed from it is trimaximal. It has four complex parameters, which can be used to 
fit the three complex neutrino masses and the remaining mixing parameter. In particular, 
if fir = fii" = fi2i then the mass matrix is also 2-3 symmetric [9], and the corresponding 
mixing is bimaximal. The mass matrix in that case coincides with the mass matrix of S4, 
with the 2 IR of S4 decomposing into (1', 1") in At, and the corresponding mixing matrix 
tri-bimaximal. In summary, the A4 effective Hamiltonian H contains one more complex 
parameter than S4, and that can be used to accommodate a non- vanishing reactor angle #13 
and a CP phase. If that parameter is chosen so that the reactor angle vanishes, then the A4 
symmetry is accidentally enlarged to an S4 symmetry. 
The neutrino mass matrix for S 3 is 

/ Hi H 2 (J>2 \ 



M„ 



(6) 



H2 A*2 A*i 

\/U 2 fii fi2 J 

where fii is the Yukawa coupling of 1 x 1 to A = 1, and fii is the Yukawa coupling of 
2 x 2 to A = 1 divided by y/2. Similarly, fi2 and J12 are the coupling constants for A = 2 
coming from 1x2 and 2x2 respectively. This mass matrix is 2-3 symmetric, befitting a 
bimaximal mixing appropriate to S 3 . It has four adjustable complex parameters that can be 
used to fit the three complex neutrino masses and the remaining complex mixing parameter. 
In particular, if fj, ± + ii 2 — fii + fi2, the mass matrix becomes magic, the mixing becomes 
tri-bimaximal, and attains the enlarged accidental symmetry S4. 

Finally let us consider S4. It differs from S4 in that (^) for 3' does not have to vanish 
(see the last row of Table 1). However, since the 3' x 3' — > 3' coupling is antisymmetric on 
the left, it cannot contribute to M u , so the mass matrices in S4 are identical to the mass 
matrices in S4. In other words, it automatically predicts a zero reactor angle 6^3, which is 
why only 3 complex parameters are enough to fit all the low energy data. 
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By using less parameters than those listed in the last column of Table 2, one can construct 
models with 'predictions' relating some of the low-energy quantities, e.g., mixing angles in 
terms of mass ratios. It should also be noted that local dynamical models may contain more 
parameters than those listed in Table 2, but if that is the case, some of them must combine 
in their corresponding LH effective Hamiltonian. 

In conclusion, if tri-bimaximal neutrino mixing is exactly true, then S4 and S4 are the 
most economical symmetry group because it makes the most predications and requires the 
smallest number of parameters to fit the remaining data. On the other hand, if 6> 13 7^ 0, 
then only A 4 among these four could be the correct symmetry group, assuming of course 
that trimaximal mixing still remains valid with better data. If neither trimaximal nor 
bimaximal mixing is valid in the face of better data, then these four groups can only lead to 
an approximate fitting in tree order. Which of the four gives a better approximation then 
depends on the data and the results of higher-order calculations. 

V. QUARK MIXING 

Quark mixing is very different from neutrino mixing because it is small. In principle 
the method reviewed in Sec. 2 can be used to find its horizontal symmetry group Q, but in 
practice, it is difficult to make Q a finite group because there is no known parameterization 
in this sector equivalent to the tri-bimaximal matrix. For Q to be a finite group, there must 
be an integer n so that g n — 1 for every g e G. If this property is true for some matrix 
g, it can no longer be true by suitably altering the matrix elements a little bit within their 
experimental error. This is similar to saying that 22/7 is a rational number but 7r is not, 
although the former is an approximation to the latter. If a finite group exists, it is likely 
that this group belongs to a series with an adjustable parameter so that the small Cabibbo 
angle can be accommodated. These would suggest finite groups like the dihedral groups D m , 
or the SU(3) finite subgroups A(3m 2 ) and A(6m 2 ). Indeed D 7 and D u have been proposed 
as candidates [1, 10]. 

However, if we want quarks and leptons to be united, like in a grand unified theory, 
then they should have the same horizontal symmetry group Q. There are two ways to 
accommodate that. One is to find a large group containing both, such as SU(3), which 
at an early stage breaks down to the individual symmetry groups in the quark and the 
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lepton sectors. Such a construction would have to be completely dynamical and the present 
approach can say very little about it. The other is to assume the leptonic symmetry group, 
like £4, A±, or S3, to be also the symmetry group of the quarks, with the same set of Higgs 
and the same vacuum alignments. Many dynamical models have been constructed along 
these lines. In that case, quarks do not mix in the tree level, so all their mixing must come 
from loop corrections. One way to do that is to assign both the LH and the RH fermions to 
the defining 3-dimensional representation, and introduce two sets of Higgs, and ip, so that 
(f> have the same vacuum alignments as $, and ip has the same vacuum alignments as ^ . The 
Higgs <fi would couple the LH fermions to the RH fermions to produce the diagonal Dirac 
mass matrices Mj (i = u, d, e, is), and ip would couple only to the heavy Majorana neutrinos 
to produce the mass matrix M N . This would guarantee Mj = M/M; to be diagonal, and 
the neutrino effective mass matrix M v = Mj M^ 1 M„ to give rise to a neutrino mixing 
appropriate to Q. How realistic schemes like that are depend on how successful are the high- 
order corrections. These corrections can be implemented by higher-dimensional terms, at 
the expense of introducing new parameters, but it would be nice if they could be computed 
in a renormalizable theory. 

I thank Christoph Luhn for helpful discussions and comments. 
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